Asymmetric dark matter may alter the evolution of low-mass stars and brown dwarfs 
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We study energy transport by asymmetric dark matter in the interiors of very low-mass stars and 
brown dwarfs. Our motivation is to explore astrophysical signatures of asymmetric dark matter, 
which otherwise may not be amenable to conventional indirect dark matter searches. In viable 
models, the additional cooling of very low-mass stellar cores can alter stellar properties. Asymmetric 
dark matter with mass 4 < Mx/GeV < 10 and a spin-dependent (spin-independent) cross section 



of a- 



10" 



10" 



cm ) can increase the minimum mass of main sequence hydrogen 



burning, partly determining whether or not the object is a star at all. Similar dark matter candidates 
reduce the luminosities of low-mass stars and accelerate the cooling of brown dwarfs. Such light dark 
matter is of particular interest given results from the DAMA, CoGeNT, and CRESST dark matter 
searches. We discuss possibilities for observing dark matter effects in stars in the solar neighborhood, 
globular clusters, and, of particular promise, local dwarf galaxies, among other environments, as well 
as exploiting these effects to constrain dark matter properties. 
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Overwhelming evidence indicates that a form of non- 
baryonic matter constitutes the majority of mass in the 
Universe. The unknown nature of the dark matter (DM) 
is a fundamental problem in cosmology and particle 
physics. Among DM candidates, weakly-interacting mas- 
sive particles (WIMPs), particularly the lightest super- 
partners in supersymmetric theories, have garnered the 
most attention [ij. We study the effects of asymmetric 
dark matter (ADM) on the evolution of very low-mass 
(VLM) stars and brown dwarfs (BD). We show that stars 
with masses Af, < O.ISM© and BDs just below the min- 
imum mass for a hydrogen-burning main sequence (MS) 
star, M > 0.05Mq, may have their evolution significantly 
altered by the accumulation of DM. 

ADM has a relic asymmetry, so it does not annihilate 
as thermal relic WIMPs do. Consequently, ADM is not 
amenable to indirect detection via observations of annihi- 
lation products from astrophysical sources |^ . Our result 
suggests a future indirect identification method for ADM 
in astrophysical sources. In the present context, the relic 
asymmetry allows stars to collect large amounts of DM 
without the accumulation being moderated by annihila- 
tion 0-01 ■ ADM models offer a possible explanation for 
the DM and baryon densities of the Universe being of 
the same order |6l4la|. ADM models are relevant to DM 
particles with masses < 15 GeV, lower than typical 
WIMPs, and interest in ADM has been fueled by possible 
direct detection signals indicating low-mass DM by the 
DAMA [13, CoGeNT [H [l^ and CRESST-II [2^ ex- 
periments |2ll . [2^ (though challenged by Xenon-10 limits, 
Ref. [H). 

It has long been recognized that DM could accumulate 
in the Sun (and other stars) subtly alteringitsproper- 
ties, particularly solar neutrino fluxes 

BSiSil- The 

scenario is simple. As the star orbits in the DM halo. 



some DM particles scatter off stellar nuclei and become 
bound to the star. The captured DM particles can be 
non-negligible contributors to energy transport in the 
stellar interior. 

DM with Afx ^ 15 GeV is captured by stars at a rate 



c: 



SI,SD 



-0.4 GeV cm 

2 



r) (l 



618 kms 



10- 

/270 kms"^ 



where M* is stellar mass, px is the DM density in 
the star's vicinity, Cp is the cross-section for ADM- 
proton scattering, Wcsc is the stellar escape speed, and 
V is the typical speed at infinity of infalling DM parti- 
cles. 
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The coefficients w 7 x 10^2 s" 
10^1 (5GeV/Afx)s"i give the rates in the Sun due 



and C|D 



to spin-independent (SI) and spin-dependent (SD) scat- 
tering respectively [25;]. We include the Aix depen- 
dence in Cq^. Cq^ is a weak function of M^. In the 
SI case, we take the scattering cross section on a nu- 
cleus with mass and mass number A to be ctn = 



A- 



where nip is the 



proton mass. We provide Eq. ([!]) for convenience, but 
calculate time-dependent capture rates using the full for- 
mulae from ^25] as described in Ref. [28| . 

VLM stars and BDs are interesting DM laboratories for 
several reasons. First, dark matter capture rates depend 
upon stellar structure only insomuch as Vcsc does. In 
VLM stars, mass is proportional to radius so VLM stars 
have similar escape speeds to the Sun and capture nearly 
the same amount of DM per unit mass. Second, stellar 
luminosity scales with mass roughly as L oc for VLM 
stars, so a Af* O.lAf© star radiates only L ~ 10~^Lq, 
and DM needs to transport a relatively small energy flux 
to alter the stellar evolution. Third, VLM stars have core 
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temperatures many times lower than the Sun and nuclear 
burning rates in this regime are rapidly- varying functions 
of temperature, so small temperature changes have large 
effects on luminosity. These considerations suggest that 
DM energy transport may dramatically alter VLM stars. 

The mean free path of captured DM exceeds stel- 
lar radii, so DM energy transport is non local. Nev- 
ertheless, we can make an order-of-magnitude estimate 
of the energy flux transmitted by DM from the stel- 
lar core by considering this transport to be diffusive, 
with an effective diffusion coefficient r/ ^ ?7((Tp/crc)^, 
where rj ~ ( 1 / rzp gp ) \/kTc/ is the standard for dif- 
fusive transport [29[. The cross-section for which the 
star is optically thick is CTc = (wp/M,)7ri?J. The factors 
(cp/cc)^ account for the facts that DM orbits span far 
less than a mean free path and DM particles scatter only 
after many orbits. Using 77' and estimating the temper- 
ature gradient as dT/dr ^ Tc/i?*, the rate at which DM 
removes energy from the stellar core is. 
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Eq. ([2]) includes only SD scattering for simplicity, iV^ 
is the number of captured DM particles, and A^p is the 
number of protons. A luminosity of L* ~ 10~'^Lq is 
typical for the « O.IM© star we have assumed. 

The energy flux carried by DM is of order the stel- 
lar luminosity (~ 10" 
cross section of cr^° = 



'Lq) if ^ 10- "iV, 



10^3. A 



10 
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cm yields a capture rate of 
Cx ^ 5 X 10^^ s"^. Assuming a lifetime r ^ 10^" yr, then 
Nx ^ 10*'^, so large effects are possible in ADM mod- 
els. We have computed the DM heat transfer rate within 
n = 3/2 polytropic stellar models, appropriate for fully 
convective stars with ~ O.IM© [30|, and confirmed 
this approximate result. These approximations are not 
self-consistent; however, they provide strong indications 
that ADM effects may be non-negligible. DM effects are 
dramatic in environments with large capture rates. We 
parametrize environment with Fb, the ratio of the cap- 
ture rate of DM to the standard capture rate in the solar 
neighborhood. As we discuss below, environments with 
Fb S> 1 are known to exist and Fb ^ 1 may effectively 
be realized in the solar neighborhood. 

We perform a self-consistent stellar evolution calcu- 
lation by computing capture rates as in Ref. 



281 and 



including DM heat transport in the Modules for Experi- 
ments in Stellar Astrophysics (MESA) software [sij. We 
base our models on the MESA VLM star models (§ 7.1 of 



Ref. [3l|). Ideally, one would calculate DM energy trans- 



port by solving the Boltzmann equation at each stage of 
stellar evolution [s^l, but this is computationally chal- 
lenging and beyond our present scope. We adopt the 



approximations of Ref. [24| , namely, that DM has an ef- 



fective temperature Tx, and that energy transport can be 
estimated from the first moment of the Boltzmann equa- 
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FIG. 1: The evolution of central temperature Tc as a function 
of central density pc for VLM stars with masses labeled at left. 
Solid lines are standard stellar models. Dashed lines include 
dark matter cooling. The line at which temperature is four 
times the Fermi temperature is marked, below which electron 
degeneracy pressure dominates pressure support. Triangles 
mark ages of 10^ yr and 5 x 10^ yr. Trajectories with one 
triangle exhibit no evolution during this period and are on 
the stellar MS. 

tion. In the case of SD scattering off of hydrogen, the 
energy per unit mass transmitted from baryons to DM is 



SD 



TT (Mx ■ 



Toe 



mpTx 



(T - Tx), 



(3) 

where is the local DM number density, /h = 0.71 is 
the local mass fraction in hydrogen, and is fixed by 
requiring that there be no net energy transfer. In the 
case of SI scattering, one sums Eq. ^ over all nuclei. 
Eq. [3] has known shortcomings, but corrections are gen- 
erally of order unity and there is no general treatment 
that can be incorporated into a stellar evolution model 
in a computationally- feasible manner Uncertainty 
associated with DM parameters is large, so this approx- 
imation should suffice for the present purpose. 

In the interest of brevity, we quote detailed results 
for a fiducial model of ADM with = 5 GeV and 

^SD 



10- 



well below current limits 331. How 



ever, we obtain similar results for SI cross-sections a^^ ~ 
10^"*° cm^ due to the large cross sections for DM scatter- 
ing on heavier elements, primarily He, C, O, Ne. Specif- 
ically, a model with Mx « 7 GeV and af 10-*° cm^ 
yields stellar evolution similar to our fiducial case and is 
of interest given indications of low-mass DM scattering 
by the DAMA [17] and CoGeNT collaborations. 
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We account for possible loss of DM by evaporation from 
the stellar interior following Ref . [s^l , with the result that 
evaporation is negligible for DM masses > 3 GeV. 
This is slightly smaller than the evaporation mass from 
the Sun (Mx ^ 3.7 GeV), primarily due to the cooler 
interiors of VLM stars. 

Fig. [T] shows the evolution of central density pc, and 
central temperature Tc, for VLM stars and BDs. This 
classic plot illustrates that a collapsing gas cloud of suffi- 
ciently low mass will achieve a maximum too meager 
to ignite hydrogen burning at a level that can halt grav- 
itational collapse j35| . Maximum Tc is achieved when 
pressure support becomes dominated by electron degen- 
eracy. Objects with > O.OSMq halt contraction when 
hydrogen burning begins and enter the stellar main se- 
quence, enjoying long MS lifetimes. Lower-mass objects 
continue to contract and cool, becoming BDs. In Fig. [1] 
we compare standard stellar evolution to evolution in- 
cluding DM cooling in our fiducial ADM model with a 
capture rate boosted by Fb = lO'^. In the DM models, 
the additional cooling causes degeneracy to he achieved 
at lower densities and temperature maxima are reduced 
at fixed stellar mass. Consequently, the minimum mass 
for MS H-burning increases by ^ 15%. 

We emphasize, with regard to Fig.[l] that the primary 
effect of ADM cooling is to drive the stellar core to de- 
generacy at a lower density. Once degeneracy sets in, 
DM energy transport becomes significantly less impor- 
tant and the degenerate core cools and contracts as in 
the standard evolution of a degenerate object [s^, [s^ . 
At all stages of evolution the stellar temperature pro- 
file is monotonic and the DM temperature closely tracks 
the baryonic temperature of the stellar core. Moreover, 
convective energy transport (see 13, 3^) remains an im- 
portant channel for energy transport in the bulk of the 
stellar interiors in all of the models we have considered. 

Important consequences of DM cooling are more ap- 
parent in Fig. [21 which depicts the evolution of luminos- 
ity for collapsing objects with 0.06 < M^/Mq < O.IIMq. 
In standard models, objects with > O.OSA/q en- 
ter a long-lived phase of constant luminosity supported 
by core H-burning, the MS. Lower-mass objects cool 
and dim incessantly. In our DM model, objects with 
^ O.IMq dim continually and exhibit no constant- 
luminosity phase, so objects with 0.08 < M/Mq < 0.10 
are not MS stars. DM cooling effects are most dramatic 
just below the minimum MS mass, M* « O.OSM©. In 
the standard case, such objects have their dimming de- 
layed by non-negligible nuclear burning, yet nuclear re- 
action rates at such low temperatures are sensitive func- 
tions of temperature [sO], so DM cooling quells this H- 
burning and drives these objects to lower luminosities. 
At M = O.OSMq the object with DM cooling is a factor 
of 15 less luminous than its standard counterpart af- 
ter 10^° yr. Notice that DM has little effect on objects 
with M < O.O5M0 because these objects are never sig- 
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FIG. 2: Luminosity evolution of low-mass stars and brown 
dwarfs. Solid and Dashed are as in Fig. [T] Dotted lines show 
evolution for Fb = 10^, 10^ for the M = O.OSMq models. 
The masses of the objects are M = O.llM©, O.IM©, O.OgM©, 
O.O85A/0, O.O8M0, O.O75M0, O.OTMq, O.O6M0, O.OSM© from 
top to bottom (every other line is labeled at the upper left). 
We show power laws for cooling of a fully convective object 
at constant effective temperature (L oc t~^^''^) and a fully 
degenerate object (L oc f~^/*) for reference (e.g., Ref. [3^). 

nificantly affected by nuclear burning. Fig. [2] shows two 
additional cases of an M = O.OSA/q object with Fb = 10 
and Fb = lO'^, illustrating that DM cooling can have 
appreciable effects with more modest DM capture rates. 

Figure [3] displays the shifts in the positions of stars 
on the HR diagram at several values of stellar mass af- 
ter 10 Gyr of stellar evolution. The changes in effective 
temperature and luminosity are significant for a range 
of stellar masses from 0.05 — 0.10 M0. ADM accelerates 
the cooling and dimming of these stars, giving rise to a 
dearth of relatively luminous VLM stars and BDs. 

It is possible to approximate the deficit of VLM stars. 
Assuming a fixed initial mass function and stellar metal- 
licity, the change in the luminosity function of VLM stars 
and BDs is determined by the mass-luminosity relation, 
L{M). We have approximated L{M) using stellar mod- 
els computed for masses separed by AM = 0.005 M© 
between M = 0.04Mq and M = 0.12Mq. Figure H 
shows the ratio of the luminosity function in our fidu- 
cial ADM scenario to a standard luminosity function for 
a stellar population after 10 Gyr of evolution. The sig- 
nificant deficit of VLM stars and BDs with luminosities 
2 X 10^^ L0 < L ^ 10~^Lq is apparent. Furthermore, 
an excess number of dim BDs near L ^ 10^^ L© is also 
evident. This pile-up results from the accelerated dim- 
ming and cooling low-mass objects in the ADM model 
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FIG. 3: Shifts in stellar positions on an H-R diagram for our 
fiducial dark matter model. Squares are standard evolution 
while triangles are the fiducial ADM model. We show shifts 
evaluated at six masses to avoid clutter in the diagram. Lines 
connect the standard and ADM cases at each mass. 

that would otherwise be more luminous VLM stars and 
BDs. 

The effects of DM on VLM stars and BDs may have 
interesting consequences. To be sure, observing VLM 
stars and BDs is challenging; however, several operat- 
ing and forthcoming astronomical facilities count obser- 
vations of local and distant VLM stars and BDs among 
their science drivers including PanSTARRS 
[381], Euclid [391], TMT |47|], GMT i48|J, and JWST 



37|, LSST 
0. 

At minimum, we have demonstrated that DM may have 
non-negligible effects on VLM stars and BDs, rendering 
them cooler and dimmer than otherwise at fixed age and 
chemical composition. The luminosity function of VLM 
stars and BDs should become significantly shallower from 
10"'^ ^ L/Lq > 10~^ and steepen at lower luminosities 
as indicated by Fig. [2] and Fig. |4l It may also be possi- 
ble to identify seemingly anomalous VLM stellar or BD 
companions in studies of transits of more luminous stars. 
The strength of this effect should be correlated with en- 
vironment: it is dramatic in regions of high-DM density 
and small in regions of low-DM density. This correlation 
may distinguish DM influences on stellar evolution from 
uncertainties in stellar modeling. 

Dramatic effects generally require capture rates larger 
than standard rates in the solar neighborhood. High cap- 
ture rates are possible in the solar neighborhood in par- 
ticular models. A co-rotating disk of DM near the Galac- 
tic plane is a prediction of hierarchical galaxy formation 



tion can enhance capture rates equivalent to Fb ~ 10^ 
or greater in the Sun and nearby stars [^, [1^. This 
suggests that observations of local VLM stars and BDs 
may constrain models of self-interacting ADM, though 
such constraints require careful modeling of stellar atmo- 
spheres and an exploration of stellar parameters, such 
as metallicity. Additionally, trace populations of VLM 
stars that were either (1) liberated in the disruption of 
merging dwarf galaxies during hierarchical formation of 
the Milky Way or (2) members of the halo star popula- 
tion that orbit through the Galactic Center would still 
be effected by ADM accumulated while in these higher- 
density environments and be subject to the DM effects 
we describe. 

Globular clusters (GC) are interesting environments 
in which to seek DM effects. The formation, evolution, 
and DM content of GCs are still under debate. However, 
GCs are many orders of magnitude denser than the so- 
lar neighborhood and have internal velocity dispersions 
a factor of 30 less than the local value [i^. More- 
over, the bounds on their DM contents are weak, with 
the ratio of DM-to-stellar mass Mum/Ms < 1 'ilj. Con- 
sequently, values of Mum/M^s well below those within 
reach of stellar kinematical studies would lead to values 
of Fb cx Px/v 1 [see Eq. ([T])]. This suggests that 
DM-induced alterations to stellar evolution may occur in 
GCs. If consensus ever formed around evidence for low- 
mass ADM, observations of stellar populations in GCs 
may reveal their past or present DM contents. 

For fixed DM parameters, high capture rates should be 
realized in particular environments. One such environ- 
ment is a Milky Way dwarf satellite galaxy. Consider the 
closest dwarf satellite. Segue I. Stellar kinematics indi- 
cate that the DM density in Segue I is 10 — 10^ times 
the local value, while typical velocities are 10^ times lower 
[ist . suggesting Fb oc p^/v 10^ - 10"'. VLM stars in 
dwarf galaxies are senstive to low-mass ADM and may 
be within reach of future observatories such as JWST. 

Interestingly, it may be possible to identify DM effects 
in the spectra of specific galaxies or GCs. Ref. 46[ em- 



and may provide a boost of Fb ^ 10 [4l|, |42|. Alterna- 
tively, DM with a significant self-interaction cross sec- 



phasized that spectral features due to Na and FeH are 
produced only by VLM stars and that the prevalence of 
these features reveals the contribution of VLM stars to 
the total galactic light. This suggests that stellar popu- 
lation synthesis studies may bear on the identification of 
DM and vice versa. Particular objects with kinematical 
or other evidence for high Fb could be targeted for obser- 
vation and the relative light contributed by VLM stars 
could limit the influence of DM on VLM stellar evolution. 
Alternatively, though more speculatively, it is thought 
that the first generations of stars formed in the inner re- 
gions of early-forming dark matter halos, environments 
with significantly higher DM density than should accom- 
pany contemporary star formation in the Milky Way, so 
galaxies in which a significant amount of star formation 
occured at high redshift may have had low-mass stellar 
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FIG. 4: Relative stellar/BD luminosity function for a stellar 
population at an age of 10 Gyr for our fiducial ADM model, 
^ai3m{L), compared to the luminosity function in the case 
of standard stellar evolution, $std {L) , assuming an identical 
IMF. 

populations markedly altered by DM. 

Several advances must be made before DM constraints 
based on this effect could be realized. First and fore- 
most, observations of VLM stars and BDs must be 
improved. Theoretically, it is necessary to implement 
compiitationally-efficient models of DM energy trans- 
port, to account for uncertainties in VLM stellar atmo- 
spheres, and explore stellar parameters, such as metallic- 
ity. We will report on these efforts in a detailed follow- 
up study. That ADM may have potentially-observable 
effects on VLM stars and BDs opens up intriguing new 
avenues for constraining DM and learning about the en- 
vironments of stellar populations. 
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